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Kondo physics in nonequilibrium interacting nanoscale devices is an attractive fundamental many- 
particle phenomenon with a rich potential for applications. Due to enormous complexity its clear 
and flexible theory is still highly desirable. We develop a physically transparent analytical theory 
capable to correctly describe the Kondo effect in strongly interacting systems at temperatures close 
to and above the Kondo temperature. We derive a nonequilibrium Keldysh field theory valid for 
a system with any finite electron-electron interaction which is much stronger than the coupling of 
the system to contacts. Finite electron-electron interactions are treated involving as many slave- 
boson degrees of freedom as one needs for a concrete many-body system. In a small vicinity of the 
zero slave-bosonic field configuration weak slave-bosonic oscillations, induced by the dot-contacts 
tunneling, are described by an effective Keldysh action quadratic in the slave-bosonic fields. For 
clarity the theory is presented for the single impurity Anderson model but the construction of the 
Keldysh field integral is universal and applicable to systems with more complex many-body spectra. 

PACS numbers: 72.15.Qm, 73.63.-b, 72.10.Fk 



I. INTRODUCTION 

With the emergence of technological feasibility in man- 
ufacturing small artificial systems with few electrons and 
using single molecules as electronic devices, quantum 
nonequilibrium has taken a new wind lasting already for 
more than two decades. On one side the diminutiveness 
of these nanoscale systems facilitates their fast passage 
to a nonequilibrium state already at tiny external volt- 
ages which calls for more comprehensive nonequilibrium 
theories beyond linear response. On the other side at 
these scales electron-electron correlations start playing a 
significant role leading to phenomena which are funda- 
mentally absent in the noninteracting counterparts. A 
remarkable representative of such no nequ ilibrium inter- 
acting phenomena is the Kondo effectP^] revealing itself 
in the differential conductance as a zero- bias anomaly 
having a universal temperature dependence^' with the 
energy scale known as the Kondo temperature, Tk. 

The Kondo physics turns out to be highly complex. 
In particular, in its temperature dependence it has 
two regimes, weak and strong coupling ones, with the 
crossover taking place in the vicinity of the Kondo tem- 
perature. This has given rise to the development of two 
types of analytical theories valid either deep below or 
much above Tk- In the construction of such theories the 
concept of a slave-boson has played an extremely fruit- 
ful role and it has been used both for the firstP^ and 
seconcP types of analytical theories. The Kondo effect 
is governed by infrared physics. A powerful method cap- 
turing this behavior is the renormalization group method 
applied to the Kondo model, both in its analyticaP and 
numerical implem entat ions. The analytical RG theory 
has been generalize d 11 ' 12 ' to nonequilibrium and provides 
a controlled theory at bias voltages eV > kT^- 

Numerical approaches have certain advantages. For 
example, the numerical renormalization group theory^ 
can cover both the weak and strong coupling regimes. 



Unfortunately, its generalization to nonequilibrium is 
still problematic. Another nu meric al method, the non- 
crossing approximation (NCAjSESl, is able to deal with 
nonequilibrium situations but it gets less reliable below 
Tk- However, the main disadvantage of all numerical 
techniques is that they do not provide a clear picture of 
a physical phenomenon and, therefore, analytical theo- 
ries have higher fundamental priority. 

The aim of the present study is to demonstrate that 
the field-theoretic approach^ based on the Keldysh field 
integral in the slave-bosonic representation is a power- 
ful tool in the development of analytical theories for the 
Kondo effect in nanoscale systems and that the deriva- 
tion of the effective Keldysh action represents a univer- 
sal, straightforward and physically clear way to possi- 
ble theoretical generalizations inspired by modern ex- 
periments. In addition to the technical flexibility this 
field-theoretic framework in its original simplest formu- 
lation is already close to the crossover between the weak 
and strong coupling regimes: it is valid for temperatures 
T > Tk- Therefore, after a proper generalization it may 
penetrate the strong coupling regime and become an an- 
alytical theory of both the first and second type. 

To be specific, here we will develop a generalization to a 
particularly complicated and at the same time extremely 
important aspect related to a strong but finite electron- 
electron interaction in a quantum dot (QD). Up to now 
many analytical theories (which are mainly of the first 
type) have only been able to describe the Kondo physics 
in QDs with either rather wealP^ or infinitely strong^HH 
electron-electron interactions. Thus, an analytical theory 
which could be placed between these two extreme cases is 
in great demand, especially, if this theory is of the second 
type which is more interesting for applications. It is our 
goal here to present such a theory. 

The paper is organized as follows. In Section [TT] we for- 
mulate the problem and transform it into a slave-bosonic 
representation. The effective Keldysh action is derived 
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in Section III It is used to derive the tunneling density 



of states in Section ITVl With Section M we conclude. 



II. HAMILTONIAN AND SLAVE-BOSONS 

We will show a nonperturbative, both in the QD- 
contacts coupling and electron-electron interaction, field- 
theoretic construction for the case of the single impurity 
Anderson model (SIAM) of an interacting QD, 
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where e d is the single-particle electronic (a =t,i) en- 
ergy level of the QD, U is the strength of the electron- 
electron interaction and n dia = d\d„ with cfi a , d a being 
the QD electronic creation and annihilation operators. 
The Hamiltonian given by Eq. ([!]) is sufficient to show 
the basic steps of the theory without technical compli- 
cations playing no principal role in the formalism. The 
generalization to more complicated models is straightfor- 
ward as will become clear from the derivation below. 

To complete the formulation of the problem one intro- 
duces the tunneling Hamiltonian, 
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dtT* 



and the contacts Hamiltonian, 
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In Eqs. ^ and ^ c\ (c a ) are the contacts electronic 
creation (annihilation) operators with respect to the state 
a, which is a quantum electronic state in the left (L) or 
the right (R) contact. The quantity T a(7 in Eq. ^ is the 
probability amplitude for the tunneling event from the 
QD state a to the contact state a, and e Q in Eq. ^ is the 
contact electronic single-particle energy corresponding to 
the state a. 

Before constructing the Keldysh field theory we rewrite 
the QD Hamiltonian, Eq. ([!]), and the tunneling Hamil- 
tonian, Eq. ([2]), using slave-bosons and new fermions re- 
placing the original electrons of the QD. To do this we will 
employ the approach developed for magnetic impurities 
by Barneses It was also used later for superconducting 
compounds by Zou and Anderson^. This approach rep- 
resents a natural extension of the infinite-?/ slave-boson 
approaches!. An alternative finite-[/ slave-boson method 
was proposed by Kotliar and Ruckensteiri 2 -^. The essen- 
tial difference between the two methods is that the lat- 
ter one, instead of new fermions, uses magnetic bosons. 
However, as mentioned in Ref. 1211 it does not matter 
how one introduces slave-particles if the corresponding 
constraints are taken into account exactly: the final re- 
sults must be independent of the slave-particle scheme. 

To apply the Barnes and Zou and Anderson technique 
one exploits the knowledge of the many-particle spectrum 



of the SIAM Hamiltonian, Eq. Q: 

£f = e l = £ d; I tl>, e n = 2e d -+ 



|0), e = 0; |t>, II), 
U. Introducing the 
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slave-bosonic creation and annihilation operators e', e 
for the empty (|0)), d\ d for the doubly occupied (| 
QD states (e-state and d-state for later reference) and 
new fermionic creation and annihilation operators p\ , p a 
for the singly occupied (| |), 
dl = epl+adtp-a, 
can be rewritten as 



= e^pa 



D) QD states, so that 
-crdpl CT , Eqs. (]J and ^ 



Ha = 



J2^plPa + C2€ d + U)dU, 



H T = Y,[Taacl(^p a + adpl 

a a 

+ T* a {epi + <rSp^ a )c a }. 
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The full Hamiltonian is H = H d + Hq + Ht- Since the 
physical subspace, the QD Fock space, is generated by 
the four states, |0), | t)i | l)i | t4)j the total number of the 
slave-bosons and new fermions, Q = e + ^2 a pl-p^ + d) d, 
must be equal to one, Q = I. One takes this constraint 
into account by replacing the total Hamiltonian with a 
different one, H — > H + fiQ, and taking the limit /1 — > 00 
at the end of all calculations. In comparison with mean 
field theories this projection onto the physical subspace 
is exact. Therefore, it is also applicable for high temper- 
atures and is suitable for a development of both the first 
and second type analytical theories. 

At this point we would like to mention that, as is obvi- 
ous from above, the applicability of the Barnes and Zou 
and Anderson approach is not restricted only by SIAM, 
Eq. ([lj. It is much more universal and applicable to 
an arbitrary many-particle system for which its many- 
body spectrum has been found. Then one can, in full 
analogy with above, introduce slave-bosons for the states 
with the even particle numbers and new fermions for the 
states with the odd particle numbers. 



III. EFFECTIVE KELDYSH ACTION 

Using Eqs. ([3), Q and ^ we develop a slave-bosonic 
Keldysh field theory. To this end we construct the 
Keldysh field integraP^. The field theory is obtained per- 
forming the functional integration over all the fermionic 
fields and eventually is formulated through the effective 
Keldysh action. In comparison with Ref. 16 this action 
depends now not only on the e-state slave-bosonic field, 
which we denote as x(i), but also on the d-state slave- 
bosonic field, which we denote as £(i): 



+ 4°iK cl (t),e q (i)]+5 t u„[x cl (t),x q W;r 1 w,e q w], 
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where x° 1< ' q ' ) W, £ cl ^ q H^) are the classical (quantum) com- 
ponents of the slave-bosonic fields, the actions and 
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Sg°j describe the free dynamics of the e- and d-states, 
respectively, and the action Stun embodies the complex 
tunneling dynamics of both the e- and d-states. 

Any QD observable O — F{d\, d a ) admits the Keldysh 
field integral representation in terms of the effective 
Keldysh action Eq. 

(6)(t) = lim e^x 

x /l?[ X (i),CW]e* Srff[xCl(t) ^ q(t); « Cl(t) '« q(t)1 x ( 7 ) 



where AFq is the normalization factor (see Ref. [TO]) and f3 
is the inverse temperature. 

The general structure of Stun is the same as in Ref. 
US S tun [x cl (t),xHt)^ cl (t),^(t)} = -ifttrln[/ + rG<°>]. 
However, the matrices and T are now different. The 
free Green's function matrix G^ has now an explicit spin 
dependence, 
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where G^ c are the fermionic Keldysh Green's function 
matrices for the isolated QD and contacts. These ma- 
trices are diagonal in the spin space. In Eq. Q we 
have picked out the spin index from the state a, i.e., 
a = {a, er} and a is the state a but the spin. This 
explicit spin dependence appears because the mapping 
— > {e^, e; d\ has an explicit spin de- 

pendence. What is more crucial is that now, due to this 
spin-dependent mapping, the tunneling matrix T is not 
anymore diagonal in the spin space. It has the same block 
structure as in Ref. [11] but its blocks have a different ar- 
chitecture in the spin space: 



M T (aat\a't') = -5(t-t r ) 
n 
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where T Serer /(t) = T a<7 a'xHt), T 5ffff /(i) = -cr'T So ._ CT /|(t). 
Here T So - CT < = 5 aa >T which assumes a symmetric coupling 
to the contacts. The matrix x(t) is the same as in Ref. 
1161 and the new matrix £ (i) is 



£(*) 
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One can see from Eqs. |9]) and ( |To| that the off- 
diagonal blocks in the spin space have an anomalous ma- 
trix structure in the Keldysh space: the classical com- 
ponents are the off-diagonal elements while the quantum 
components are on the main diagonal. Physically this 
non-diagonal spin structure and the anomalous struc- 
ture in the Keldysh space are the consequence of the QD 
population fluctuations between one and two electrons. 
These fluctuations are now allowed due to the finite value 
U of the electron-electron interaction. 



At this point it is important to mention that in our 
field theory there is no need to scan over a huge fam- 
ily of diagrams in order to find diagrams recovering the 
symmetry between virtual transitions into the empty and 
doubly occupied states. Violation of this symmetry rep- 
resents the main problem in extensions^ of NCA for 
finite-l/ systems. In contrast, in our field theory on one 
side one does not have to battle with the diagrammatic 
combinatorics and on the other side one does not need 
to care about the symmetry violation since both virtual 
transitions are equally treated by keeping in the effective 
Keldysh action the same powers of both the empty and 
double occupancy slave-bosonic fields. Below we keep 
terms only up to the second order in these fields and it 
is obvious that the same can be done at any order. From 
the diagrammatic perspective the theory could be close 
in spirit to the equilibrium imaginary-time approach^ 
for the infinite-?/ Anderson model. 

As we are going to construct an analytical theory of 
the second type, we consider a QD in the Kondo regime 
formed by weak (see Ref. ffij]) slave-bosonic oscillations 
excited by the electronic tunneling between the QD and 
contacts. In this case one can expand Stun around the 
zero configurations of the slave-bosonic fields %(t) and 
£(f) up to the second order: 



Stun[x c \t) lX q (t);e\t),e(t)] 

_ 7lT 
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dt' (e l (t) p(t)) x 



where E^ ,K and E^ ,K are the slave-bosonic self- 
energies of the e- and d-states. The self-energies Egg _,K 
have been derived in Ref. Q1)] and the new objects, 



X 
X 
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where x = L, R and the functions gt' :K , g£, g ± are the 
same as in Ref. [16l We use the Lorentzian density of 
states for the contacts, ^c( e ) = v c W 2 /(e 2 + W 2 ), and 
define T = 4ttv c \T\ 2 . 

The advantage of this quadratic model is that it is 
analytically solvable. However, its applicability is limited 
by temperatures T > Xk and QD chemical potentials 
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T < fio — ^ U — r which assumes that the theory is 
valid for U ^> T. The latter is exactly our primary goal. 

Let us say a few words about the structure of the tun- 
neling action. As one can see from Eq. (JTlJ, it does not 
contain terms mixing \(t) and In the second order 
such terms just do not appear. The physical picture be- 
hind this is that the charge excitations corresponding to 
the single and double occupancies do not interact. This 
is valid for large U. However, even for large U the ex- 
citations can interact at low temperatures through the 
Kondo resonance for fj, in the vicinity of the symmetric 
point, e d + U/2. Therefore, one expects that in this vicin- 
ity the quadratic theory is not valid at low temperatures. 



IV. TUNNELING DENSITY OF STATES AND 
RESULTS 

It is now a straightforward task to calculate observ- 
ables using Eq. ( |lT| ). The observable storing the whole 
universe of the Kondo physics in QDs is the tunneling 
density of states (TDOS), u a (e) = -(l/S7r)Im[G| (rCT (e)], 
where G d aa , (e) is the retarded QD Green's function. Us- 
ing the functional integral representation, Eq. ([7|, and 
performing the Gaussian integral we obtain 



v a U) = Z(e) 



[e d -e + ff ER( e )P + [ ff Si(e) 



r + 



[e d + U-e + gY%{e)] 2 + [gY} d (e)] 2 



(13) 



where g = nT/2. In Eq. (Il3j Se (I) is the real 



(imaginary) part of the e-state retarded slave-bosonic 

self-energy (see Ref. [16| and £^ is the real (imag- 
inary) part of the d-state retarded slave-bosonic self- 
energy, £*(e) = We/[n(e 2 + W 2 )] - E?(e), E d (e) = 
-W 2 [2 ~ n L (e) - n R (e)]/[2Tr(e 2 + W 2 )], where n L , n R are 
the Fermi-Dirac distributions of the contacts electrons, 
n L(R)( e ) = { ex P[/?( e — Mo ± e ^/2)] + l} -1 , characterized 



by an external voltage V. The function Z in Eq. (13 1 is 



Z(e) 



x 



r w 2 

An e 2 + W 2 

»(e d )[l - n B d(e d )} 

1 - "L(ed)«R(ed) - n Bd (e d )n(e d ) ' 



(14) 



where n(e) = n L (e) + n R (e) - 2n L (e)n R (e), n Bd (e) = 
{exp[/3(2e + f/-2 Mo )] + l}- 1 . 

The QD TDOS, Eq. Q, is shown in Fig. [l] Since 
U is finite, the TDOS has two charge excitations, the 
single and double occupancy excitations close to e d and 
e d + U, respectively. For a QD in which the strength U 
of the electron-electron interaction is finite the behavior 
of the Kondo resonance as a function of /i — £d is quite 
different from the limiting case U = oo. For a finite 
value of U the Kondo peak first decreases when (j,q — 
e d increases up to the symmetric point U/2 and after 
this point increases again when /io — e d approaches U, 
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FIG. 1: (Color online) The analytical result, Eq. (131, for 
the QD TDOS. Here kT = 0.0035r, e d = 0, ^ - e d = 2.4r 
and 97.6r, U = WOT, W = 1000r. For nonequilibrium eV = 

o.4r. 




FIG. 2: (Color online) The universal temperature depen- 
dence of the conductance. Insets show the two qualitatively 
different regimes in the Kondo temperature dependence on 



so that two symmetrically located Kondo peaks, close to 
ed and e d + U, have the same height as demonstrated in 
Fig. [T] This means that the Kondo temperature Tk as a 
function of /Un — e d decreases up to U/2 and then increases 
again. Our Keldysh field theory correctly predicts this 
behavior of Tk obtained from the scaling in the universal 
temperature dependence of the conductance as shown in 
Fig. [2] This result agrees with the general expression 
(see, e.g., Ref. kT K /T ~ (W/T) exp{-27r0uo-ed)[#- 
(Mo — e d)]/r?7}. In Fig. [2] we also compare our theory 
with the numerical renormalization group theory^^. 
This comparison proves that, indeed, our theory is of the 
second type and it is reliable for temperatures T 2Tk- 



5 



V. CONCLUSION 

In summary, we have developed an analytical nonequi- 
librium Keldysh field theory for the Kondo effect in QDs 
with finite electron-electron interactions (U) which are 
much stronger than the QD-contacts coupling (r). The 
theory is nonperturbative in both U and T and valid for 
temperatures T ^ 2Xk- Although, for clarity, it has been 



presented for the case of SIAM with normal contacts, the 
construction of the Keldysh field integral has a universal 
theoretic scheme applicable to setups with ferromagnetic 
or superconducting contacts coupled to any interacting 
nanoscopic system after its many-particle spectrum has 
been found. 
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